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ANNOTATED  COMPUTER  OUPUT  FOR  SPLIT  PLOT  DESIGN:  SAS  GLM 


W.T.  Federer,  Z.D.  Feng,  M.P.  Meredith,  and  N.J.  Miles-McDermott 
Mathematical  Sciences  Institute,  Cornell  University,  Ithaca,  NY 

ABSTRACT 

fr' 

The  analysis  of  covariance  for  split  plot  designs  is  not  always 
straightforward  when  using  a  statistical  software  package  such  as 
SAS  PROC  GLM.  In  order  to  demonstrate  correct  analyses  several  data 
sets  are  examined  and  annotated  SAS  output  is  given.  Hypothetical 
data  are  analyzed  first  without  and  then  with  the  covariate  included. 
The  whole  plots  are  arranged  in  a  RCBD  and  the  covariate  is  measured 
on  the  subplot  experimental  units.  A  second  example  has  whole  plots 
arranged  in  a  CRD  and  the  covariate  measured  only  on  the  whole  plot 
experimental  units. 

Complete  ANOVA  tables  for  both  examples  may  be  computed  in  a 
single  procedural  call  to  SAS  PROC  GLM.  Both  Type  I  and  Type  III 
sums  of  squares  are  necessary  to  construct  the  proper  ANOVA  table. 

A  commonly  employed  approach  requiring  two  separate  procedural  calls 
to  GLM  is  also  demonstrated.  Formulae  for  the  standard  errors  of 
the  difference  between  adjusted  whole  plot  and  subplot  means  are 
reported . 


INTRODUCTION 

This  is  part  of  a  continuing  project  that  produces  annotated 
computer  output  for  the  analysis  of  balanced  split  plot  experiments 
with  covariates.  The  complete  project  will  involve  processing  three 


Ls» 


mare 


examples  on  SAS/GLM,  BMDP/2V,  SPSS-X/MANOVA ,  GENSTAT/ANOVA ,  and 
SYSTAT/MGLH Only  univariate  results  are  considered.  We  show  here 
the  results  from  SAS  GLM . 

For  Example  1,  the  data  are  artificial  and  were  constructed 
for  ease  of  computation;  the  experiment  design  for  the  whole  plots 
is  a .randomized  complete  block  and  the  split  plot  treatments  are 
randomly  allocated  to  the  split  plot  experimental  units  within  each 
whole  plot.  Example  2  is  the  same  as  Example  1  except  that  a 
covariate  varies  from  split  plot  to  split  plot.  The  data  for 
Example  3  come  from  an  experiment  wherein  the  whole  plot  treatments 
are  laid  out  in  a  completely  randomized  design  and  the  split  plot 
treatments  are  randomly  allotted  to  the  split  plot  experimental 
units  within  each  whole  plot.  The  value  of  the  covariate  varies 
from  whole  plot  to  whole  plot  but  is  constant  for  all  split  plots 
within  a  whole  plot  treatment. 

We  present  the  elementary  computational  steps.  Simple 
hypothetical  data  are  used  for  the  first  two  examples  so  that  it  is 
easy  to  provide  all  detailed  computations  to  illustrate  hew  each 
number  is  obtained.  Some  readers  may  wish  to  skip  the  detailed 
computations  (see  Federer,  1955,  Chapter  XVI).  The  third  example 
comes  from  Winer  (1971).  The  detailed  computations  are  given  in  his 
book  (p.  803) . 


2 


Data  S P-1 


Split  plot  data  with  whole  plots  arranged  in 
randomized  complete  block  design 
(hypothetical  data) 


Wl 

Whole 

plot  treatment 

W2 

split 

plot 

treatment 

split  plot  treatment 

Block 

S1 

S2 

S3 

S4 

Total 

S1 

S2 

S3 

S4 

Total 

1 

3 

4 

7 

6 

20 

3 

2 

1 

14 

20 

2 

6 

10 

1 

11 

28 

8 

8 

2 

18 

36 

3 

6 

10 

4 

4 

24 

10 

8 

9 

13 

40 

Total 

15 

24 

12 

21 

72 

21 

18 

12 

45 

96 

Total  and  Means 


Blocks 

(8  observations) 


_ Total _ Mean 

140  5 

2  64  8 

3  64  8 

Grand  Total  168 

Grand  Mean _ 7 


W(whole  plots) 
(12  observations) 

Total _ Heap 

Wl  72  6 

W2  96  8 


S (split  plot) 

(6  observations) 

Total _ Mean 

36  6 

42  7 

24  4 


66  11 


Model : 


Yijk 


=  u  +  p  .  +  t.  +  6.. 
3  1 


+  “k + 


<“T)ik  +  £ij* 


ji  =  mean 

Pj  =  effect  of  block 

<5 ^ j  =  error  (a) 

€ijk  =  error  <b) 


j 


(cxt) 


ik 


=  effect  of  whole  plot  i 
=  effect  of  split  plot  k 
=  effect  of  interaction  of 

whole  plot  i  and  split  plot  k 


where  it  is  assumed  that  p^  - 
and  Py  6y,  and  are 
j=l,2,***,r,  and  k=l,2,***,s. 


N(0'CTp )  '  >  N<o.«rj>, 

mutually  independent. 


ftijk  "  N(°'cre) 


i=l,2, 


,a 


3 


(*)Notation  follows  that  of  Searle (1971) ;  since  the  design  is  balanced, 
R(p \n,T,a,aT)  -  R ( p | 4 ) ,  etc.  The  simpler  notation  is  used  later. 

(**)  SxB:W  means  SxB  within  W. 

Calculations  of  SS's: 

N  =  2*3*4  =24  ,  Y  =  7 

2  3  4 

R(M,p,T,6,a,ar(e)=  212  Y?  . .  =  (32+  62  +  62+***+  182  +  132)  =  1616 

i=l  j=i  k=l 

R(p)  =  NY2  =  24 • (7) 2  =  1176 
R(p,T(6,a,aTfe |p)  =  1616  -  1176  =  440 

2  2  2 

R  (p  |p)  =  R  (p ,  p )  -  R(p)  =  +6g  +—  1  -  1176  =  1224  -  1176  =  48 

2  2 

R(-rlfi)  =  R  (p ,  t  )  -  R(p)  =  ^72~^96  J  -  1176  =  1200  -  1176  =  24 
R(6|p,p,T)  =  R(6,p,p,T)  -  R(p,p)  -  R(t,p)  +  R(p) 


2  2  2  2  2  2 
[20  +28  +24  +2  04  +  3  6<5  +  40‘5' 


-  1224  -  1200  +  1176 


=  1264  -  1224  -  1200  +  1176  =  16 

R (a  |p)  =  R(a,p)  -  R(p)  =  ,(31  +6§  )  _  1176  =  1332  _  1176  =  156 

R(ar  |jj,a,r)  =  R(ar,p,a,r)  -  R(p,a)  -  R(m,t)  +  R(p) 


22222222 
15  +24^+12  +21  +21  +18  +12  +45^' 


-  1332  -  1200  +  1176 


Nl 


B1 

1 

2 

1 

2 

6 

B2 

2 

2 

0 

4 

8 

B3 

3 

5 

2 

0 

10 

Total 

6 

9 

3 

6 

24 

=  1440  -  1332  -  1200  +  1176  =  84 

R(£  \n,P  ,6  ,a,r  ,ar)  =  R  (e  ,\x  ,a  ,  p  ,  6  ,  t  ,  ar  )  -  R  (u,p,t,6)-  R(fi  ,a  ,t  ,ar)  +  R(r,p) 

=  1616  -  1264  ~  1440  +  1200  =  112 

Data  SP-2 

Data  SP-2:  Data  SP-1  with  the  following  covariate  Z 
which  varies  with  split  plot 

Covariate  (Z) 


whole  plot 
W1 _ 


Total 


S3  S4 


Total 


9  15 


blocks 

(8  observations) 


W  (whole  plot) 
(12  observations) 


S  (split  plot) 
(6  observations) 


■ijLLLeM 

Me 

1 

14 

14 

2 

20 

20 

3 

26 

26 

Model :  Y .  . ,  =  (j+p  ,+r  .+o  .  ,+a.  +  (ar)  . . +0 .  (Z  .  .  -Z  )  +p  _  ( Z .  • ,  -Z  .  .  ) +e  .  . , 

ijk  H  l  lj  k  '  ' lk  Klv  ij*  •••'  2 '  ijk  l j • '  xjk 

=  whole  plot  regression  slope  P2  =  split  plot  regression  slope 

where  p,  Py  <5^,  ak,  (aT)ik-  and  are  as  in  SP-1,  and  Z^.  and 

are  the  arithmetic  means  for  Z..,  . 

•  •  •  ijk 


Table  of  sum  of  squares  and  cross  products 


Source _ df _ XI _ YZ _ £Z 


B 

2 

48 

18 

9 

W 

1 

24 

12 

6 

BxW  (error  a) 

2 

16 

4 

1 

S 

3 

156 

33 

9 

SxW 

3 

84 

33 

21 

SxB:W  (error  b) 

12 

112 

17 

20 

Mean 

1 

1176 

420 

150 

Total 

_ 2A _ 

1616 _ 

537 

216 

YY  column  is  the  same  as  in  SP-1,  ZZ  column  is  computed  in  the  same 
fashion.  Thus,  only  computations  for  YZ  column  are  illustrated. 


2  3  4 

Total  =2  2  2  Y..  -Z... 

YZ  i=l  j=l  k=l  13k  1Jk 

=  3(1)+  6(2)+***+  14(4)+  18(4)+  13(7)  =  537 

Meanyz  =  NY...Z...  =  =  420 


3  2  4  2  4 

2(2  2  Yiik)(2  2  Z..  ) 

~  1=1  1=1  k=l  J  i=l  k=l 

Byz  -J - 


2*4 


-420  -40(14)+64(20)+64(26) 
8 

=  438  -  420  =  18 


-  420 


2  3  4 


3  4 


.2  (.2  2  Y..  ) (  2  2  Z..  ) 

W  =  .21  ^  _  420  = 


432  -  420  =  12 


3(4) 


2  3  4  4 

2  2(2  Y  .  . .  )  (  2  Z  .  ..  ) 

i=l  -i  =1  k=l  k=l 

BxWyz  =  ■  J  1  K  X - — -  438  -  432  +  420 


=  454  -  438  -  432  +  420  =  4 


2  3 


2  3 


4  ‘A  A  Yijk}  (  A  AZijk> 


s  =  z  1=1  3  =  1  Zl  1=1  ^  =  1 

°Y7  ^ 

k=l  2(3) 


-  420  =  453  -  420  =  33 


2  4  3  3 

A  2  <  2  YijA  2  zijk) 

SxWyz=  -  - J— i - 3— = -  453  -  432  +  420 


=  498  ~  453  -  432  +  420  =  33 

6 


SxB:W  =  2  12  Y.  ..Z...  “  454  “  498  +  432 

YZ  i=l  j=l  k=l  1Dk  13k 

=  537  -  454  -  498  +  432  =  17 


Source 


B  (block) 

W  (whole  plot  treatment) 
Regression  (a) 

BxW  (error  (a) ) 

S  (split  plot  treatment) 

Sxw  (interaction  of  S  and  W) 
Regression  (b) 

SxB:  W  (error  (b) ) 


Total  (corrected  for  mean) 


=  R(p |4,t) 

=  R (t | 4,p,P1) 

=  R(P1Ip,p/t) 

=  R(6  Im,P,t,P1) 

=  R (a  \u  ,P  ,t  ,clt  ,P2) 

=  R(ar \u,p, t ,a,P2) 

=  R(P2ln»P/T»a>aT) 

=  R(t  \  u  ,  p  ,a  ,t  ,clt  ,P_) 


3.4286 

16.0 


84.243 

37.474 

14.450 

97.550 


P2  =  SxB:Wyz  /  SxB:Wzz  =  17/20  =  0.85 

The  SS's  adjusted  by  regression  on  Z  are  illustrated  below: 

R(p|u)  =  48,  remains  same  since  it  is  not  of  interest  to  adjust  blocks 
for  Z. 


R(t,5|»x,p,P1)  =  (Wyy+  BxWyy)  - 


(Wyz+  BxWyz) 
WZZ+  BxWzz 


=  (24  +  16)  -  2  —  =  40  -  =  3.4286 

6  +  1  7 


R(6  |4,p,T,P1)  =  BxW^ 


i^aL  .  -  a!  .  o 

BxW„„  1 


R(t  |»if  P  =  R(r,<5  UifPfP^  -  R(6  |p  ,  p  ,r  ,PX) 


2  56 

=  40  -  -  0  =  3.4286 


«*t.r 


■> 


.1] 


(Bxw  )-  2 

R(P,  |u,T  ,p)  =  -  =  —  =16 

B*WZZ  1 


R(tt/£.|fl/P,T,CtT,P_)  (  Syy  +  Sx  B  !  Wyy  )  ~ 


(Sy„  +  SxB:W  ) 


SZZ  +  SxB:Wzz 


=  (156  +  112)  - 


9+20 


=  268  -  86.207  =  181.793 


R(aT,e \n,p,a,T,P  )  =  (SxW  +  SxB:Wvv)  - 


(SxWv„  +  SxB:Wv  ) 


SxWZz  +  SxB:Wzz 


=  84  +  112  -  x  '  =  196  -  60.976  =  135.024 

21+20 

Note:  R  (a  ,  e  |  |i  ,  /3  2 )  and  R(ot  ,e  |p  ,a  ,  r  ,/?2 )  are  intermediate  steps  for 

later  use. 

(SxB:W  ) 2  2 

R(P,  |p  ,P  ,<*  ,+  ,ar)  =  -  =  -  =  14.450 

z  SxB:Wzz  20 

(SxB:W  )2  172 

R  ve  Im  ,  p  ,a  ,t  ,ar  ,0  )  =  SxB:W  -  -  =112  -  -  -112  -  14.45  -  97.55 

z  XI  SxB:Wzz  20 

R  (a  |  p  ,  p  ,  t  ,aT  ,  0  )  =  R(cr,t  |fi,p,T,ar,0  )  -  SS  error  b  =  181.793  -  97.55 

-  84.243 

R(ar  |p,p,a,r,P2)  =  R  (<xt  ,  e  |  p  ,  p  ,a  ,t  ,  02 )  -  R  (e  |  p  ,  p  ,a  ,  r  ,ar  ,  02 ) 

=  135.024  -  97.55  =  37.474 


isavyvv>aviviv»v.vivw.<.l,i 


•7*>  '>  ">  V  V 


Data  SP-3 

Split  plot  data  with  plots  arranged  in  a  completely 
randomized  design  and  a  covariate  Z  that  is  constant 
within  the  whole  plot.  (Winer,  1971,  p.  803) 


whole  plot  Subject 


Split  plots 


Total 

Y 


5 

6 

7 

8 

Total 

Mean 


10 

8 

3 

18 

15 

12 

5 

27 

20 

14 

8 

34 

12 

6  . 

2 

18 

15 

10 

1 

25 

25 

20 

8 

45 

20 

15 

10 

35 

15 

10  . 

2 

25 

132 

95 

39 

227 

16.5 

11.9 

4.88 

Model : 


yijk  *  “  +  Ti +  6ij  +  ak  +  <Ta)ik +  Pi<zij-z..>  +  njk 
feet  (whole  plot)  6..  =  error  (a)  e. ..  =  error  (b) 

1 J  1 J  K 

=  whole  plot  regression  slope 


t ^  =  A  effect  (whole  plot) 
=  B  effect  (split  plot) 


where  6^  -  N(o,o|),  ~  N(o,c*),  6^  and  are  mutually 

independent.  i=l,2,**‘,a,  j=l,2,***,r,  and  k=l,2,***,s. 


Analysis  of  variance  and  covariance 
Source _ _ 


A  (whole  plot)  =  R(t|h,0  ) 

Regression  =  R(0  |m,t) 

Error  (a)  =  Rfftlu,!-,/^) 

B  (split  plot)  =  R(a|pi,T,ar) 

AxB  (interaction)  =  R (ra | 4 , t , a ) 


Error  (b) 


=  R  (t  |m,t  ,a,Ta) 


Total  (corrected)  =  R(t  ,a  ,p,  ,tcx  ,  6  I4) 


SS 

44 .492 

166.577 

61.298 

85.563 

0.563 

6.375 

388.438 


*  ** 


I*  I*1 


Control  language  Is  typed  in  upper  case  and  comments  are  bolded. 


DATA  ONE; 

INPUT  BLOCK  WHOLE  SUBPLOT  Y;  =>  Input  variables 

TITLE  SP-1 :  SPLIT  PLOTS  WITH  WHOLE  PLOTS  ARRANGED  IN  RCB  DESIGN? 

CARDS ;  =>  Tells  SAS  that  data  follow 

1113 

112  4 

113  7 

:  =>  Data  are  entered  with  only  one  datum  per  line 


3  2  4  13 
PROC  GLM ; 

CLASS  BLOCK  WHOLE  SUBPLOT;  =>  Designates  classification  variables 
MODEL  YIELD=BLOCK  WHOLE  BLOCK*WHOLE 

WHOLE  WHOLE* SUBPLOT/ SS 3  P;  Designates  model  being  used.  The 

SS3  option  requests  only  type  III  sums 
of  squares  and  P  requests  residuals 
(only  one  type  SS's  was  requested 
because  the  data  are  balanced  making 
all  types  SS's  equal).  Type  I  SS's  are 
the  cheapest  to  compute. 

TEST  H=BLOCK  WHOLE  E=BLOCK*WHOLE ;  =»  Requests  SAS  to  test  the  whole 

plot  effects  using  error (a) 

Note:  SAS  always  computes  F  tests  based  on  the  residual  sum  of 

squares.  This  is  not  always  the  appropriate  test  in  split  plot 
analyses  so  adding  the  TEST  statement  (above)  is  critical  to  obtaining 
an  appropriate  test. 


Note:  Because  estimates  of  both  a  whole  plot  regression  slope  and 

split  plot  regression  slope  are  needed,  two  procedural  calls  to  SAS  GUM 
are  required.  The  first  call  gives  the  appropriate  whole  plot  analysis 
and  the  second  gives  the  appropriate  split  plot  analysis. 

Procedural  Call  for  Whole  Plot  Analysis 

DATA  ONE; 

TITLE1  SP-2 :  SPLIT  PLOT  DESIGN  WITH  WHOLE  PLOTS  ARRANGED  IN  RCB: 

TITLE2  WITH  A  COVARIATE  VARYING  WITH  SPLIT  PLOT; 

INPUT  BLOCK  WHOLE  Z1  Z2  Z3  Z4  Y1  Y2  Y3  Y4  ; 

Z= (SUM (OF  Z1-Z4 ) / 2 ) ;  =>  Z  and  Y  are  scaled  for  this  analysis 


Y= (SUM (OF  Y1-Y4 ) /2 ) ; 
CARDS ; 


1 

2 

3 

1 

2 

3 


4  7  6 
10  1  11 
10  4  4 
2  1  14 
8  8  2  18 
10  8  9  13 


PROC  GLM; 


TITLE3  CORRECT  WHOLE  PLOT  ANALYSIS; 

CLASS  BLOCK  WHOLE; 

MODEL  Y=BLOCK  WHOLE  Z/SOLUTION  SSI  SS3  P; 


so  that  the  sums  of  squares  are  correct 


For  whole  plot  analysis  data  must  be 
organized  in  a  similar  arrangement 
with  all  split  plot  values  for  a 
particular  BLOCK  by  WHOLE  combination 
on  the  same  line  (see  INPUT  statement 
above  for  order) 


LSMEANS  BLOCK  WHOLE/STDERR; 
ESTIMATE  'WHOLE  PLOT  SLOPE'  Zl; 


The  SOLUTION  option  yields 
the  parameter  estimates  and 
so  gives  the  estimated 
regression  slope 

Yields  adjusted  treatment  means  and 
standard  errors. 

The  ESTIMATE  statement  gives  the 
estimated  regression  slope  and  its 
standard  error  directly. 

Procedural  Call  for  Split  Plot  Analysis 

DATA  TWO; 

INPUT  BLOCK  WHOLE  SUBPLOT  Z  Y; 

CARDS ; 

11113 

112  2  4 

113  17 


3  2  4  7  13 
PROC  GLM; 

TITLE3  'CORRECT  SPLIT  PLOT  ANALYSIS'; 
CLASS  BLOCK  WHOLE  SUBPLOT; 

MODEL  Y=BLOCK  WHOLE  BLOCK* WHOLE  SUBPLOT 
WHOLE*SUBPLOT  Z/SOLUTION  SSI  SS3  P; 
LSMEANS  SUBPLOT  WHOLE*SUBPLOT ; 

ESTIMATE  'SUBPLOT  SLOPE'  Z  1; 


Note:  Even  though  we  are  estimating  only  one  slope  in  this  example, 

two  procedural  calls  are  required  in  order  to  estimate  the 
regression  slope.  If  both  the  whole  plot  and  split  plot  are 
specified  in  one  run,  Z  becomes  confounded  in  SUB (A)  and  01 

cannot  be  estimated  correctly. 


Ii  W  ttBH  MS  BBi 


DATA  ONE? 

INPUT  A  Z  Y1  Y2; 

SUBJECT  =_N_; 

MY= (SUM (OF  Y1-Y2) )/(SQRT(2) ) ;  =>  Z  and  Y  rescaled  so  SS's  agree  with 

those  of  Winer 

Z  =  2*Z/ ( SQRT ( 2 ) ) ? 

TITLE1  SP-4 :  SPLIT  PLOT  DESIGN  WITH  WHOLE  PLOTS  ARRANGED  IN  CRD; 
TITLE2  WITH  A  COVARIATE  CONSTANT  IN  SPLIT  PLOT? 

CARDS ; 

1  3  10  8 
1  5  15  12 
1  8  20  14 


2  2  15  10 
PROC  GLM ; 

CLASS  A? 

MODEL  MY=Z  A/SOLUTION  SSI  SS3  P? 
TITLE3  CORRECT  WHOLE  PLOT  ANALYSIS; 
LSMEANS  A/STDERR; 

ESTIMATE  ' REGR  SLOPE'  Z  1; 


rocedural  Call  for  Correct  Split  Plo 


DATA  TWO; 

INPUT  SUB  A  B  Y; 
CARDS ; 

1  1  1  10 
112  8 
2  1  1  15 


8  2  2  10 
PROC  GLM; 

CLASS  SUB  A  B; 

MODEL  Y=A  SUBJECT (A)  B  B*A/SS1  SS3  P 
TITLE3  CORRECT  SPLIT  PLOT  ANALYSIS; 
LSMEANS  B  B*A/STDER R; 


Variances  and  Standard  Errors  of  Adjusted  Means  and  Differences 

Amongst  Adjusted  Means  for  SP-2 


/ar<Yi.k  adj>  =  (ap  +  °l  +  +  K  +  s*6><Zi..-Z...>2/W>‘Bzz 

+  al  (zi.k'zi..)2/SxB:Wzz 

/ar(Yi--  ad  j  ^  =  [al  +  s(CTp  +  CTfi)3/rs  +  (a*  +  Z...)2/WxB2 

^ar(Y-*k  ad  j }  =  (ap  +  a6  +  °2)/ar  +  al  (z..k‘  Z...)VSxB:Wzz 

( Z  .  ,  -  Z  .  )  2 

,ar<*i--  adj  -  V--  adj>  -  -  “6>  fi  +  1  WXB,/" 


,  „»  I 


l/ar(Y  ,  ,  .  -  Y  ,  ,  , . )  =  a  —  + 

'  **k  adj  ‘-k  ad j '  e  ar 


(z..k--  z..k) 

SxB: W„„ 


"ar<Yi.k  adj  -  Yi-k'  adj  ^  = 


2  2  ^Zi*k‘”  Zi*k^ 

°a  r  +  SxB:  W,,- 


and,  for  i  *  i ' 


Var ( Y .  .  .  .  -  Y . ,  .  ,  . . )  =  -  (a2  +  a? )  + 

v  i*k  ad]  l  *k  ad] '  r  '  e  6' 


(zr..-  Z...) 

WxB_„ 


(at  +  Sa5) 


(Z..  ,  ,  -  Z  .  ,  -Z.  .  +  Z  .  ) 

'  i  • k  i  • •  i*k  i»  • ' 


SxB:  Wr 


Estimates  of  the  variance  components  or2  and  ct2  are  required  to 

calculate  standard  errors  of  the  above  differences  amongst  adjusted 
treatment  means.  From  the  expected  mean  squares  of  the  ANOVA  table  it 

is  known  that  error  (a)  and  error  (b)  estimate  o2  +  s a2  and  a2  , 

respectively.  If  error(a)  and  error(b)  are  denoted  Efl  and  E^, 

respectively,  then  a2  is  estimated  by  (Eg-  Efa)/s  •  Hence,  the  desired 
standard  errors  are  given  by: 


•  i*  •  aaj  1  •  *  aaj  •  J  3 


SE^Y**k  adj  ~  Y**k'  adj)  I  Eb  ar  +  SxB:W( 


L  1 2  .  (Z..k'“  Z*.k} 

—  I E,  I _ +  7Z77+rm  —  ~ 


„  2  ,  (Zi*k'“  Zi*k^3 

—  E.  —  +  - „  ■  - 


SE<Yi*k  adj  -  Yi.k'  adj)  -  Eb  ?  +  ±SxB?W 


and,  for  i  *  i’ 

,2[E  +(s-l)E  ]  (Zi...-  Z...)2 

SE(Yi-k  adj  "  Yi-k-  adj)  "  \-^Ti - -  +  —ttr -  Ec 


(Zi'*k<_  Zi‘-.“  Zi-k+  zi--)2  V  ^ 2 

+  SxB:WZZ  Efc^ 


Variances  and  standard  errors  of  adjusted  means  and  differences 
amongst  adjusted  means  for  SP-3. 


Var  (Y .  . 

1  i*  k 

adj ) 

-  K  +  °l 

Va*-(Yi.. 

adj ) 

=  ( CT 2  +  SO 

£ 

Var(Y..k 

adj ) 

-  K  +  °l 

Var  (Y . 

x*  • 

adj 

-  Y . 

l  • •  adj 

Var (Y  , 

adj 

Y,*k’  adj 

V3r(Yi.k 

adj 

Yi*k'  adj 

(Z.  -  Z  ) 

'  5  •  •  •  * 


E  (a) 


(Zi.~  z..) 

E(a)„ 


<Zi'.-  Zi>) 

E (3) Tf 


and  for  1  *  i'  , 


Var<*i-k  adj  -  adj>  =  2K  +  °l^r  +  K  +  "V 


zz 


Estimates  of  the  variance  components  o*  and  o*  are  given  by  E(b)  and 

[E(a)  -  E(b)]/s,  respectively.  Hence,  the  desired  standard  errors  are 
given  by: 


SE<*i..  adj  '  *i'..  adj) 


E  (a) 


(Zi*.’  Zi-> 


sr 


E  (a) 


ZZ 


\ 


SE(Y--k  adj  *  Y..kadj>  =  J^(b)/ar  =  SE(Y..k-  Y..k.) 
SE<5i-k  adj  -  5-.k’  adj»  *  J^ETb)??  "  SE(Y1.k-  Yi.k.) 


and,  for  i  *  i'  , 


SE^i-k  adj  -  V.k‘>  = 


2  fE  fa)  +  ( s-1)  E  fb) 1  (Zi'»  Zi- } 

E  (a) 


sr 


2  \ 

-  Efa) 


ZZ 
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